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We present an extension of relativistic single-particle distribution function for weakly interacting 
particles at local thermodynamical equilibrium including spin degrees of freedom, for massive spin 
1/2 particles. We infer, on the basis of the global equilibrium case, that at local thermodynamical 
equilibrium particles acquire a net polarization proportional to the vorticity of the inverse temper- 
ature four-vector field. The obtained formula for polarization also implies that a steady gradient of 
temperature entails a polarization orthogonal to particle momentum. The single-particle distribu- 
tion function in momentum space extends the so-called Cooper-Frye formula to particles with spin 
1/2 and allows to predict their polarization in relativistic heavy ion collisions at the freeze-out. 

£ . I. INTRODUCTION 

' The single-particle distribution function is the main quantity in kinetic theory and its form at local thermodynamical 
, equilibrium for relativistic, weakly interacting, gases is well known. For spinless particles, it is simply the Bose-Einstein 
y—i • distribution with x-depcndcnt values of temperature and chemical potential, also known as Jiittner distribution: 

co ; 

co i f( x >p) = e f}( X y p -z( X ) _ i M 

co : 

f^) , where p — -^rU is the inverse temperature four-vector, To being the proper temperature measured by a comoving 
thermometer with the four- velocity u and £ = fio/To is the ratio between the proper chemical potential fig and To- 
The above formula has, as straightforward consequence, the invariant momentum spectrum at local thermodynamical 



equilibrium, the so-called Cooper-Frye formula 
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where E is a space-like 3-dimensional hypersurface. This formula is widespreadly used in e.g. relativistic heavy ion 
collisions to calculate hadronic spectra at the end of the hydrodynamical stage. 

The Jiittner distribution ([1]), multiplied by a degeneracy factor (25 + 1), is also used for particles with spin (for 
fermions replacing the —1 with +1 in the denominator) being understood that / means the total particle density in 
phase space, i.e. summed over polarization states. However, in general, particles may not evenly populate the various 
polarization states and one may then wonder what is the appropriate extension of ([I]) in this case. Indeed, in this 
work, we will answer this question and provide a generalization of ([1} and ^ including the spin degrees of freedom. 
We will argue that a non-even population of the polarization states arises when the inverse temperature four-vector 
field has a non- vanishing antisymmetric part of its gradient and calculate the polarization vector for massive spin 1/2 
particles. Phenomenologically, this extension may have several interesting applications. For instance, it would make 
it possible to predict the value of particle polarization in relativistic heavy ion collisions 043 at the hydrodynamical 
decoupling, provided that local thermodynamical equilibrium applies to spin degrees of freedom as well. 

The derivation of the extension of Cooper-Frye formula will be done in several steps and it also requires a summa- 
rization of relativistic kinetic theory from quantum field viewpoint including spinorial degrees of freedom. The reader 
who is only interested in the final result may jump to sect. I VIII 
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The paper is organized as follows: in sect. [TT] we set the stage for the generalization of the single-particle distribution 
function according to relativistic kinetic theory in a quantum field framework; in sect. IIIII we summarize results 
obtained for global thermodynamical equilibrium with rotation deriving more compact formulae for spin 1/2 particles; 
in sect. IIVI we generalize the single-particle distribution function for spin 1/2 particles to local thermodynamical 
equilibrium and discuss the physical meaning of the antisymmetric tensor coupled to the spin matrices; in sect. IVll we 
calculate the polarization vector in a relativistic fluid and in sect. lVIll we summarize the results and draw conclusions. 



Notation 



In this paper we adopt the natural units, with h = c = K = 1. 
The Minkowskian metric tensor is diag(l, —1, —1, —1); for the Levi-Civita symbol we use the convention e 0123 = 1. 
We will use the relativistic notation with repeated indices assumed to be saturated. Operators in Hilbert space will 
be denoted by an upper hat, e.g. R, with the exception of the Dirac field operator which is denoted with a capital ^. 
Repeated greek indices are saturated, repeated latin indices are not. 



II. RELATIVISTIC KINETIC THEORY AND QUANTUM FIELDS 

Classically, hydrodynamics emerges as an effective description of the underlying microscopic dynamics, specifically 
of an underlying kinetic theory, where the substratum is a large number of interacting particles. In this approach, 
all hydrodynamical quantities and, chiefly, the stress-energy tensor in relativistic hydrodynamics, can be expressed in 
terms of the single particle distribution function f(x,p): 

T^(x)= J^p^p"f{x,p) (3) 

where f(x,p) fulfills a transport equation, notably Boltzmann equation. On the other hand, (relativistic) hydrodynam- 
ics can also be seen as the realization, on the average, of conservation equations holding at the more fundamental level 
of quantum field operators. In this approach, the stress-energy tensor is seen as the mean value of the corresponding 
quantum one, which is a function of the fundamental fields: 

T^ix) = tr{p:T^(x) :) (4) 

being p is the density operator (either pure or mixed) and the normal ordering is meant to remove unwanted diver- 
gencies [7|; in general, a renormalization procedure is in order when calculating mean values in an interacting field 
theory [8j. 

What is the relation between the definitions (|3l) and (|4|)? The connection between relativistic kinetic and quantum 
field theories has been nicely described in ref. [7] and the conceptual tool bridging the gap between them is the 
covariant Wigner function. For a charged scalar field "0 it is defined as: 

W(x, k) = J d 4 y 2 e-* k -y (: ft (x + y/2) $(x - y/2) :) (5) 

where (.) stands for tr(p.). Provided that the interaction terms are small and that significant W(x,p) variations 
occur on a macroscopic scale, i.e. much larger than Compton wavelength, k is an almost on-mass-shell four vector 
and the covariant Wigner function can be written: 

W+(x,k) ee 9(k°)W(x,k) = J ^ 5 4 (k-p)f(x, P ) 

W-{x,k) = e(-k°)W{x,k) = [ — S 4 (k + p)f(x, P ) (6) 



which define the distribution functions f(x,p) and f(x,p) with on-shell four- vector p and e = \/p 2 + m 2 . For a free 
field, the definition ([5]) and the ^ lead to: 
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f( X >ti = 2(2^)3 J duS ( U -P) e m ' X ( a p-u/2 a P+u/2) 
^ (X ' P) = 2(2^)3 / d4u6 ( U -P^ e ' lU ' X ( b l-u/2 b P+u/2) 
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a p , bp being destruction operators of particles with four-momentum p normalized so as to: 

[ap,aj,] =2 £ <S 3 (p-p') 

Thus: 

/ d 3 x f(x,p) = i<ata p ) = ^ jd 3 * /(*,,) = = (7) 

which is exactly what one would like to have for the non-interacting case, namely the space integral of the distribution 
functions is the number of (anti)particles per three-momentum cell. 

The function / turns out to be the familiar single-particle distribution function in phase space and, in the weakly 
interacting case, it can be shown to obey the Boltzmann equation Q. 

The covariant Wigner function makes the connection between and ^ manifest. If the quantum stress energy 
tensor of the free scalar field is chosen to be: 



then 0): 



T^{x) = (: f^{x) :} = ~(: ft (x)d^d u ^(x) :) = J d 4 fc k*WW(x, k) = J ^ p»p u [f{x,p) + f(x,p)] (8) 

where the eq. ([6]) has been used. We note in passing that a different form of the stress-energy tensor (differing thereof 
by a divergence) of the charged scalar field would have not led to the last, somewhat familiar, expression. We do not 
know whether a suitable change in the definition of W, / or both, allows to keep relations ([7} and © for any form 
of the quantum stress-energy tensor. Also, it is worth pointing out that, in general, it has been shown that different 
quantum stress-energy tensors are thermodynamically inequivalent [t| [Io| and it would not be thus surprising if a 
change of the quantum stress-energy tensor would not preserve the relation (J8j> . 

A similar connection can be built up for the spin 1/2 particles and the Dirac field. In this case, the covariant 
Wigner function is a 4 x 4 spinorial matrix: 

W(x, k) AB = - 7 ^- J f d 4 y e^ k -y(: V A (x - y/2)* B (x + y/2) :) 



(27T) 



(2tt)4 



d 4 y e- lk - y (: V B (x + y/2)# A {x - y/2) :> (9) 



VP being the Dirac field 0- We note in passing that the above definition has to be modified in full spinor electrodynamics 
d, [ll| to preserve gauge invariance, but this can be neglected for the scope of this work. 
Likewise, the distribution functions, for on-shell k, are defined by an equation similar to ([5]): 

W + (x,k) = 8(k°)W(x,k) = ± I < ^5 i {k~p)Y,Mp)frs(x,p)u s (p) (10) 

r.s 

W-{x,k) = d{-k°)W{x,k) = -\ [ ^S 4 (k+p)Y,Vs( P )frs(x,p)v r (p) (11) 



where u r (p),v r (p) are the spinors of free particles and antiparticles in the polarization state r (namely a helicity or 
third spin component) normalized so as to u r u r = 2m, v r v r = —2m. Note that the distribution functions are 2x2 
matrices, whose diagonal components f rr are the phase-space densities of particles in the polarization state r, that is 
(for the proof, see Appendix A): 

d 3 X f rr (x,p) = = ^ J d3x frr(x,p) = ^(i^A.r) = ^ (12) 



1 It should be reminded that the normal ordering for fermion fields involves a minus sign for each permutation, e.g. : aa^ 
Therefore, taking into account anticommutation relations, for fields : \Jf b(v) '■= ~ '■ ^b(v)^a( x ) '■ 
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It is now convenient to introduce a more compact notation, where / is a 2 x 2 matrix and U(p) = (u + (p),u-(p)) a 
4x2 matrix (and U(p) = C^(p)7° a 2 x 4) such that the equation (fTU|) can be written as: 

W(x,k) = J / — 5 4 (k-p)U(p)f(x,p)U(p)-S 4 (k+p)V(p)f T (x,p)V( P ) (13) 



The mean value of any operator involving a 4 x 4 matrix A such as (: \I r A3' :) can then be written, using ((9|) and (fT3 

~2s 



(: *(x)A*(x) :) = / d 4 fc tr(^(x, fc)) = / ^ [tr(tf(p)/(x,p)tf(p)i4) - tr^.F^FWA)] 



/d^p ^ 



= / ^ tr 2 (/(x,p)C7(p)A[/(p)) - tr 2 (/ T (x,p)y(p)^(p)) (14) 

We have introduced the symbol tr 2 to mean the sum over polarization states (indices r, s), whereas tr without 
subscript stands for the trace of 4 x 4 spinorial matrices. 

We are now in a position to calculate the mean value of physical densities, like currents, stress-energy tensor etc. 
The mean value of the current \D r 7^\E r , is a straightforward application of the eq. (fl"4f with ^1 = 7^: 



f{x) ee (: ^{x)j^{x) :)= —Y,frs(x,p)u{p) s 'l t "u(p) r ~f rs (x,p)v r (p)^v s (p) 



rs 



= I ^I>.(*>P)V^.-/r.(*,p)2p' i *r,= f ^p»(tT 2 f(x,p)-tT 2 f(x,p)) 

J rs J 

where we have used the well known properties of spinors u, v. The meaning of the last expression is clear: the mean 
value of the current is the momentum integral of the single-particle distribution in phase space multiplied by their 
four-velocity, with a + sign for particles and a — sign for antiparticles. In full accord with eq. (|12[) . tr 2 / and tr 2 / are 
the particle number densities in phase space. 

Similarly, one can obtain the expression of the mean value of the canonical stress-energy tensor (see Appendix B for 
the proof): 



T'-ur) ee 1<: ^(x)^d^(x) :) = J ^Vp" (tr 2 /(x,p) + tr 2 />,p)) (15) 



and, by using eq. (|14j) . the mean value of the canonical spin tensor: 

S x ^(x) ee ¥(x){ 7 \ -) = \J 5 tf2 U^P)U(P){1\ V"}U(p)) - tr 2 (f T ' (x,p)V(p){ 7 \ V"}V(p)) 

(16) 

where S^" = (i/4)[7 /i , 7"] are the generators of Lorentz transformation of spinors. Note that the mean value of the 
canonical stress-energy tensor turns out to be symmetric even though the operator is not. This implies, for the free 
Dirac field, that the divergence of the mean value of the spin tensor (which equals the antisymmetric part of the mean 
value of the stress-energy tensor because of angular momentum conservation) vanishes, that is d\ (: S <^ u :) = 0. This 
result should have been expected in a theory without interaction. We note that this equation holds only if the density 
operator is time-independent, which, as we will see, is not the case for the local equilibrium density operator. 



III. GLOBAL THERMODYNAMIC AL EQUILIBRIUM WITH ROTATION 

Atglobal thcrmodynamical equilibrium with finite angular momentum density the density operator is well known 
fl2l ITa j and reads: 

p=|cxp[-i//T + ^g/T + ^.J/T]P y (17) 

where a; is a constant fixed vector whose physical meaning is that of an angular velocity and T is the global temper- 
ature, that is the temperature of a thermostat in contact with the system or that measured by a thermometer at rest 
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with respect to the external inertial observer. The Py operator is the projector operator onto localized states [14j 
which is needed in order to avoid the relativistic singularity at r — c/cu. In this case, the single-particle distribution 
function of an ideal relativistic Boltzmann gas of particles with spin S can be obtained from just statistical mechanics 
arguments [Til ]: 



f(x,p) r 



— e^ e" 



- ' p ~ (£> s ([p]- 1 Ro,(io;/T)[p]) + D s ([p]tR (i (io;/T)[p]t-i)) i 



(18) 



where /3 = A(l, u; Xx) is the inverse temperature four- vector; A is the fugacity; [p] is the SL(2,C) matrix corresponding 
to the Lorentz transformation taking the time unit vector t into p (so-called standard transformation); D s stands 
for the (S, 0) irreducible representation of SL(2,C) (the (0, S) being D s t -1 |T5j); R is the SL(2,C) corresponding of a 
rotation, which is calculated for an imaginary angle iu/T. Note that, according to the eq. (|18[) the matrix / is not 
diagonal and: 



tr 2 s +1 /=e*e^nr 2S+1 R*(*u;/T)=e£e-^ £ e^/ T = e« ^ Pp x(^) 



We now focus our attention on the S =1/2 case. Since Weyl's representation the Dirac spinors, with the normal- 
ization u(p) r u{p) s = 2mS rs ) and with C = ia% (<x being Pauli matrices) [15l Il6j read: 



U{p) 



( D s {\p\) 
WW' 1 ] 



V(p) 



DS([p]t-i<7) 



(19) 



and D 1 / 2 (R&(ioj /T)) = exp[(— w/T)<73/2] one can rewrite the (fT8"|) . for spin 1/2 particles in the Boltzmann limit as: 

1 



f(x,p) 



2m 



being: 



S 2 = - 
2 



U(p)exp[(u/T)J: z ]U(p) 



a 3 
a 3 



whereas for antiparticles it can be shown that (see Appendix C) 

1 



-e-«e"^ 



2m 



V(p)exp[-(w/T)Z z }V(p)Y 



The eqs. f|20p and (1221) can be written in a fully covariant form by introducing the tensor: 



where tt^ turns out to be the acceleration tensor of the Frenet-Serret tetrad of the /? field lines \U 
second equality holds for a rigid velocity field only 0E3]- Hence, in the Boltzmann limit: 



(20) 
(21) 

(22) 

(23) 
The above 



f(x,p) 



e«e"^ 



-^—U{p) exp 
2m 



1 



U(p) 



f(x, P ) = -c-te-P-P 



1 

2m 



V(p) exp 



1^ 
2 



V(p) 



The extension to Fermi-Dirac statistics of this formula should, in principle, be determined calculating the covariant 
Wigner function (JSJ with the operator (fTT|) . However, this is not a straightforward calculation like in the non- rotating 
case, hence we make an ansatz about this extension which reproduces the known distribution in the non-rotating case 
and, at the same time, has the (f24)) as its Boltzmann limit. We first introduce the 4x4 matrices X, X whose relation 
with /, / reads: 



1 



U(p)X(x,p)U(p)=f(x,p) 



1 



2m ' 2m 

and assume that they are, at thermodynamical equilibrium with rotation: 

X = I exp[/3 • p — £] exp 



[V{p)X(x,p)V(p)] T = f(x,p) 



1 " 




-l 




w : n 


+ ') 






2 









X = exp[/3 ■ p + £] exp 



"1 " 




— w : n 




2 





(25) 



(26) 



6 



where : is a shorthand for the rank 2 tensor contraction, i.e. A : B = A^B^ . Consequently, the single-particle 
distribution functions read: 



f(x,p) = ^~ u (p) ( ex P[/3 ' P ~ £0*0] ex P 



1 " 




w : 2-i 


+ ') 


2 





f{x,p) = - 

which have (1241) as Boltzmann limit 



1 

2m 



V(p) ( exp[/3 -p + (,] exp 



"1 ■ 






—vj : 2j 






2 







1 T 



V( P ) 



(27) 



IV. LOCAL EQUILIBRIUM OF A RELATIVISTIC FLUID OF PARTICLES WITH SPIN 



In quantum relativistic statistical mechanics, local thcrmodynamical equilibrium density operator reads: 



1 

Ple = ^ — exp 



- / d 3 x Mx)f^(x) - \vj, v {x)S^{x) - £(z)?(z) 



(28) 



and it is obtained by maximizing entropy S = — tr[plogp] with the constraints of given local values of mean energy- 
momentum, angular momentum and charge density [18]. As entropy is not conserved in nonequilibrium situation, 
the above operator breaks covariance (there cannot be invariant spatial integrals of non-conserved currents) and it is 
time dependent. This operator is used in derivations of the relativistic Kubo formulae of transport coefficients 11911 
and, in comparison with usual formulations, it has an additional term involving the spin tensor, obtained in ref. [lOj. 
In principle, all quantities at local thermodynamical equilibrium in quantum relativistic statistical mechanics should 
be calculated using (j2"8"j) as density operator, including the covariant Wigner function of the Dirac field. However, the 
full calculation is quite complicated and goes beyond the scope of this work. At the lowest order of approximation, 
however, we know that it must yield the same formal expression at global thermodynamical equilibrium with space- 
time dependent intensive thermodynamics functions, that is space-time dependent /3, w and £. Hence, the matrices 
(f2T))) become: 



A 



exp[/3(x) • p — £(x)] exp 



1 



vo(x) : E 



X 



exp\j3(x) • p + £(x)] exp 



vj(x) : E 



(29) 



and, accordingly, the single-particle distribution functions ([271 



f(x,p) = ^U( P ) (exp[/3(z) .p-S(x)]exp[-~w(x) : E] + lj U(p) 
f(x,p) = -7^(V(P) (exp[/3(x) .p + £(z)]exp[~ W (z) :£]+/) V(p)) T 



(30) 



The matrices X and X can be expanded into a power series, each term having an even number of 7 matrices because 
E/ij/ oc [7^,71/]. Therefore, the particle phase space densities read: 



tr 2 / = —tT 2 (U(p)XU(p)) 
2m 



l-tr{XU{p)U(p)) = ±-tr(X{f + m)) = UiX 



1 



tr 2 / = -—tr 2 {V{p)XV{p)) = -— tr (XV (p)V(p)) = -—tr(X(f - m)) = -trA (31) 
where we have used trace cyclicity and the spinorial completeness relations (written in compact notation): 

U(p)U( P ) = {]/> + m)I V( P )V(p) = {f- m)I (32) 

Thus, we obtain the free Dirac part of the canonical stress-energy tensor (IT51) as: 



T^(x) 



d 3 p 



^p"(tr 2 / + tr 2 /) 



d 3 p 



p^p" [trA + trA] 



and the current: 



f{x) = \ J^P" [trA -trA] 



(33) 



(34) 
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both involving the traces of X, X. If vj, which is adimensional in natural units, is small enough such that one can 
write the expansion: 



e /3(z)- PT ^) exp 



T^w(x) : E 



j =£(-i)«+ 1 



e nf3(x)-p±n£{x) 



n=l 



±-w{x) : E 



(35) 



then it is possible to obtain an approximate expression of those traces for <C 1. Note that at full rotational 
equilibrium this condition amounts to require Tiuj/KT <C 1 (natural units purposely restored) which is a normally 
fulfilled condition. Hence: 



4 



tr I exp 



±— w(x) : E 
. 2 w 



71 77, 

tr ( I ± -ro(ir) : E + —w(x) : E ro(a;) : E ) = 4 • 



n7(x) Ap ro(x) ffT tr (E Ap E aT ) (36) 



where the tracelessness of E matrices has been used. By using known formulae for the traces of 7 matrices it can be 
shown that: 



tr (EApEo-r) — g\a9pr — 9\T9pa 



whence the eq. (|36)l becomes 



tr ( exp 



Therefore, we have: 

n=l 

where: 



±— vj(x) : E 
2 v ; 



-w(x) : zu(x) 



4 + — zu(x) : m(x) 



(37) 



4np + — rii?(l — ?t.f)(1 — 2np)tu(x) : w{x) (38) 



1 



np 



and similarly for X with the replacement £ — >■ — £. One can then plug the eq. (|38| and its corresponding for X into 
eqs. (13"3")) and to obtain, e.g., for the charge density: 



d 3 p (np — np) + vj(x) : vj(x)— I d 3 p [np(l — np)(l — 2np) — np(l — np)(l — 2np)] 



3°{x) = 

For w = one recovers the usual expression; it is worth noting that the lowest order correction to charge density and 
stress-energy tensor is quadratic in w, i.e. in huj/KT at equilibrium. 

It is now crucial to know the physical meaning of the tensor w, whose expression is, as yet, known only at equilibrium. 
While the general physical meaning of the fields j3 and £ can be easily inferred from the equilibrium limit (/? is the 
local inverse temperature four- vector field and £ the ratio between the comoving chemical potential ^to(^) an d the 
local comoving temperature Tq(x) = \/ \J]P), ro AI ^(a;)'s expression cannot be uniquely obtained from the equilibrium 
distribution. The reason of this ambiguity is that at rotational equilibrium, the tensor w is, at the same time (see 
eq. (US])): 



where f2 M „ is the acceleration tensor of the Frenet-Serret tetrad of the /3 field lines and @,[T3: 



(39) 



(40) 



In a noncquilibrium situation, the right hand sides of eqs. (|40|) and (|39|) differ and it is not obvious which one applies, 
perhaps neither. However, if the system is not too far from equilibrium, zu cannot be too distant from the right hand 
side of P0"|) . Particularly, the difference must be of the 2nd order in the gradients of the (3 field, i.e.: 



-(d^ v -d v (3 tl )+0{d 2 f3) 



(41) 



where (D(d 2 f3) stands for terms which are quadratic in the /3 gradients (for instance: (d ■ (3)(d^p v — d v (3p)) or second 
order gradients (like 0^8^/3^) which vanish at equilibrium, where [tI Il7jj: 

Pn = bp + n7 p i" 

being b and w constants. For the lowest-order formulation of relativistic hydrodynamics, the expression (|41l) is 
sufficient to determine the expression of spin-related quantities, as we will see in the next two sections. 
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V. THE SPIN TENSOR 



In order to estimate the polarization of particles, what will be done in the next section, it is necessary to calculate 
the spin tensor first. We start by working out eq. ([16]) with the help of eq. (|30|) : 



S X >^ = \ tr 2 (f{x,p)U(p)tf,'Sr'}U(p)) - tr 2 (/ T (x,p)^(p){ 7 °, E^}V(p)) 



i r d 3 P 

Am J 2e 
1 f d 3 p 

Am J 2e 
1 f d 3 p 

Am J 2e 

i r d 3 p 

A. ~2T 



tr 2 (tj{p)XU{p)U{p){ 1 \^}U( P )) +tr 2 (F(p)AV(pMp){ 7 \ TT}V(p)) 
tr 2 (C/(p)X(^ + m){ 7 \ E""}l7(p)) + tr 2 - m)X{ 7 A , E^}V(p)) 

tr + m){ 7 \ E""}^ + m)) + tr (X{f - m){ 7 A , TTty - m)) 
tr ({X,j*}{ 7 A , E M "}) - tr ({X,|*}{ 7 A , E^}) 



(42) 



where we have used the relations (|32p and taken into account that X, X are a linear combination of terms with an 
even number of 7 matrices. By taking advantage of the cyclicity of the trace, the last expression can be rewritten: 

= \ J ^ tr {{f, 7 A }{X, E^}) + tr {\ST t f\ [X, 7 A ]) - tr ({f, 7 A }{X, E^}) - tr ([£"", jf] [X, 7 A ]) 



and, using the known relations: 



i,~f x } = 2p A 



5 A '^ = \f^T Vtr ({X, E^}) - tr (i(p" 7 " - pV)[*, 7 A ]) - 2p A tr ({X, E^}) + tr (i(p^V - pV)[*, 7 A ]) 



1 fd 3 p 
~2e 



Aphi (XE^) - p^tr (i[ 7 A , 7 ,y ]X) + p K tr (i[ 7 A , ^]X) 



-4p A tr (XE"") + p"tr(z[ 7 \ 7 ,y ]X) - p"tr (z[ 7 \ 7 ^]X) 
= I / lj 4p>tr ~ 4 ^' tr ( SA " X ) + 4 ^tr (E A ^X) - 4/tr (XS"") + 4p"tr (E Aly X) - 4p"tr (E^X) 

Finally, we can write the spin tensor as: 



& 



^ _ 1 /" d 3 ? 



(43) 



where 



0"" = tr {X^ lu ) 0^ = -tr (IS"") 



Altogether, the form (|43p of the canonical spin tensor only depends on the fact that X, X are a superposition of an 
even number of 7 matrices. The full antisymmetry of the indices is now apparent, although it was already contained 
in the operator definition ensuing from the properties of 7 matrices. Thanks to trace cyclicity, can be written as a 
derivative with respect to the w tensor: 



0^ = tr(XE^) = ^(-1)™+V 



■n/3(x)-p+n£(x)± 



( 



tr I exp 



—zu(x) : E 
2 w 



= y (_ ^n+l e -n^)- P +n^) i _^L_ tr ( exp 



n=l 

d 



— mix) : E 
2 v ; 



tr log < / + c 



-P(x)- P +ax) 



exp 



-/cu(x) : E 



(44) 
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where the (135)) has been used. Then, using the approximation (|37|) . the eq. (|44]) becomes: 

G"" ~ Y(-l) n+1 e- ni:iixyp+<{x) -^—- ( 4 + — 07(2;) : rote) ) 
00 

= ^(-l) n+1 c- nl3{x) - p+n ^ x) nm{x)^ = n F (l - n F )w{x)^ (45) 



n=l 

Likewise, it can be shown that: 



~ n F (l - f^Mx)^ (46) 

Note that, according to (j4"5|). (|4"5"j) and (l4"6")l the spin tensor, at the lowest order, depends linearly on zu, unlike the 
current and the stress-energy tensor. Now, being: 



\j^T P XQ " U = J ^ P X n F (p)(l - n F (p)) = -w^r— J ^ n F (p) 

and defining the scalar function: 



e ^ 7 e e^'P-« + 1 



W^/e^-g^r (47) 



one can write: 

2 J e e ~ dp 2 ' 
and similarly for the antiparticle term. Therefore, plugging eq. (|4"T)l into (|4"5)) : 

<S A ^ = tu A n7(j:) Ml/ + rotation of indices + I u x ru(x) tiu + rotation of indices (48) 

where 

1 OF 

i = ~T W 

and similarly for I. The formula (|48p extends a formula obtained in ref. [l4j for the ideal relativistic Boltzmann gas 
at equilibrium. In fact, in that paper, it was assumed that the spin tensor S could be factorized into the product 
of the four- velocity u x and an antisymmetric rank 2 tensor a^" . We can now see that such an assumption is in 
general not valid. Even in the simplest case of an ideal Fermi gas at equilibrium the spin tensor does not factorize. 
However, the time component of the spin tensor 5 0,/il/ at equilibrium, what can be defined as "spin density", can be 
shown to agree with the expression found in ref. [l4j . Because of the axial symmetry of the operator (fTT)) . the only 
non- vanishing component of a rank 2 antisymmetric tensor, such as and is the (12) or (21), and so from (l43l) : 



5°' 12 (x) = - / d 3 p tr(X£ 12 ) - tr(X£ 12 ) (4<h 



2 . 

Taking the Boltzmann limits of X, X and rewriting £ 12 = £ 2 : 
S°' 12 {x)= l - J d 3 pe"^ {e«tr(exp[( W /r)E z ]S z )-e-«tr(exp[-( W /T)I] 2 ]E 2 )} 

= lj d 3 pe-^ | e «^tr(exp[(a;/r)Ej) + e-«^tr(e X p[-(a;/T)Sj)| 

= l (e e + e -C) J d 3 pe-^^ 5 ;tr(ex P [( W /T)Ej) = i(e«+e-f) J d 3 p e"^ ^2tr(exp[( W /T) ff 3/2]) 

= (e« + e-«) | d 3 p e-^* ^ tr(exp[(o;/r)«73/2]) (50) 

where the third equality ensues from the invariance of this trace under a change of sign of uj. This result indeed 
coincides with the spin density calculated in ref. [14] for an ideal relativistic rotating gas at equilibrium, obtained 
without the use of quantum fields. 
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VI. POLARIZATION 



The main consequence of the distribution functions (J30J) is that spin 1/2 particles get polarized at local thermody- 
namical equilibrium. In this section, we will calculate the polarization four-vector, which, for a particle with mass m 
and four-momentum p is defined as: 



1 aP T 

Iln = — —tiipcrT s p (si) 

2 m 



where S pa is the mean value of the total angular momentum operator of the single particle. If we want to know the 
mean polarization vector of a particle with momentum p around the space-time point x, we cannot but divide the 
total angular momentum density in phase space J r °< pT (x,p) by the density of particles in phase space, that is: 

/n , U 11 dJ°^( X ,p) P T 

0W*.p)> = - 2 ^j^r d3p - (52) 

where J(x) is the total angular momentum density: 

J x ' pa {x) = x p T Xa {x) - x a T XT {x)+S x ^{x) 
Replacing the stress energy tensor with its expression in (TT51) : 

, 3 [ j = W - x°p p )tr 2 f(x,p) + {X) 
a^p a^p 

The Levi-Civita tensor in the above equation makes the orbital part of the angular momentum density irrelevant, so 
we are left with the spin tensor contribution only. Thus, taking into account that the particle density in phase space 
is tr 2 / = (l/2)trX (see eq. (JSJ)) and using (|4"3")l for the spin tensor, for particles we have: 

(iW*,p)> - --^e WCTT - ( P °e p ° +p°e 0p + p p e°°) p - = ~-^—e^ T (53) 

4tr2/ e m ztiX m 



and similarly for antiparticles with replacing 0. At the lowest order in w, using eq. (|45|) and eq. (|38|) : 

(Up)(x,p) ~ ~\e, paT (1 - n F )m(x) p ^ (54) 

Finally, not far from equilibrium, using (|41[) : 

(n„(i,p)) - ^ P *T (1 - n F ) {d p P° - d°p>>) V - = \e, paT (1 - n F )d p P° V — (55) 
16 m 8 m 

We can see that for a degenerate Fermi gas, the polarization tends to zero because all levels are filled. For a non- 
degenerate Fermi gas, the formula coincides with the one obtained for the ideal relativistic Boltzmann gas 20], as it 
should. 

The polarization vector can be separated in its time and space parts as: 

n = (n°, n) = ^^((v x p) ■ P , e (v x p) - ^ x P - v/?° x P ) (56) 

The above formula has the remarkable con seq uence that quasi-free particles get polarized not only in a vorticous flow 
(what was pointed out in previous works [20(), but also in a steady temperature gradient without velocity flow, i.e. 
when V/3° ^ 0. Restoring natural constants, in the non-relativistic limit where the distinction between comoving and 
observed temperature can be neglected, the predicted polarization reads 

n = (n ,n) = (i-n F )^-^(o,VTxp) 

which is usually tiny but could be relevant in some extreme situations. 

Now, the polarization 3- vector n in the rest frame of the particle with four-momentum p can be found by Lorentz- 
boosting the above four- vector [2l| : 

n = n- g n P (57) 

eye + m) 
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implying a longitudinal polarization (helicity): 

n p = n p - p 2 n ■ p = ( x - v ) ^^e(v xp)-p = ^^(v x p) ■ p (58) 

e(e + m) \ eye + m) J 8m 8 

Hence, the mean helicity of a spin 1/2 particle is approximately proportional to the curl of the (3 = 7V/T0 field. 

It may be of interest, e.g. for relativistic heavy ion collisions, to calculate the space-integrated mean polarization 
vector. For a three-dimensional spacelike hypersurface E, one has, using (1431) : 

ds A L.(-i/ 2 ) ewwT _£- / ds A / e- 

(U ^' ]) = " = -i e -^^LY ( 59 ) 



f p 4 ?7x diV 

/ dSA — iT 2j(x,p) £ ^ 



Furthermore, taking into account the particle phase space density is tr2/ = (l/2)trX and using (|38|) at the lowest 
order in w. 



one finally obtains, using eqs. (|4"5)) and (|1D]): 



dE A p A tr 2 /~2 / dS A p np 



1 p : 



/" dE A p A n F (l - n F )n7 p<T 1 pT J dS A p A n F {\ - n F )d p ^ 



(n^(p)) ~ --e W ar JTJ ^ a e W><7T 7 ( 60 ) 

4 m diV 8 m / ,„ A 

e— y dS A p n F 

These results deserve some discussion concerned with the physical meaning of the spin tensor, which has been the 
crucial ingredient to obtain (|53l) and the ensuing formulae. In principle, the definition of the polarization vector 
(151 p involves the total angular momentum of the particle, hence the formula should be invariant under a change 
of the stress-energy tensor and spin tensor operators keeping the integral of the total angular momentum density 
invariant, a so-called pseudo-gauge transformation of the stress-energy tensor [22| (see also detailed discussion in 
refs. dElii). However, in a formula like (1551) . one defines the local value of polarization and, therefore, a dependence 
on the particular spin tensor is implied. It could be expected that a space-integrated expression like (|59]l would be 
independent of the spin tensor choice, in fact this is not the case because of the explicit time-dependence of the local 
equilibrium density operator (see discussion in sect. IIV|) . Such a dependence does not enable to write the mean value 
of the divergence of an operator as the divergence of its mean value, thus breaking the pseudo-gauge invariance of the 
total angular momentum. For instance, had we used the Bclinfante symmetrized stress-energy tensor, the ensuing 
value of polarization at local thcrmodynamical equilibrium (|59[) would vanish. To summarize, the choice of a specific 
spin tensor operator is necessary to calculate the polarization of particles and we have chosen the canonical spin 
tensor (see eq. (IT5|) . which is the same used in ref. [23| to calculate the polarization of electrons. Even though it might 
appear disturbing that polarization at local thermodynamical equilibrium depends on the particular quantum spin 
tensor (whence the stress-energy tensor) of the theory, it has been recently shown that in thermodynamics this is a 
general feature [H, [l(| ■ 

VII. SUMMARY AND CONCLUSIONS 

In summary, we have obtained the form of the single-particle invariant distribution function for particles and 
antiparticles with spin 1/2, generalizing the Jiittner distribution: 



f(x,p) rs = ^—u r {p) [ exp[/3(x) • p - £(x)] exp 
f{x,p) rs = --^-v s (p) [ exp[^(x) -p + £(x)}exp 



■ij u s(p) 
i) V r (p) (61) 



where f3 is the inverse temperature four-vector, £ is the ratio between comoving chemical potential and temperature, 
E^" = (i/4)[7 M , j v ] are the generators of Lorentz transformations of 4-components spinors, u(p) and v(p) are the 
spinors solutions of the free Dirac equation and: 

tn^ = -i(^-^/3 M )+0(a 2 /3) (62) 
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This formula leads to a generalization of the Cooper-Frye formula for particles with spin: 



(63) 




(64) 



where rip is the usual Fermi- Jiittner distribution. Therefore, particles acquire a polarization in vorticous flows, which 
has a longitudinal component and a transverse polarization even in steady temperature gradients without flow. 



This work was partly supported by HIC for FAIR and the Italian Ministery of Research and Education PRIN 2009. 
We thank L. Csernai for interesting discussions and for useful suggestions. 
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To prove eq. (fT2")l . we first manipulate (fTU|) and (fTTj) multiplying both sides by u(k) c on the left and u(k)d on the 
right: 
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APPENDIX A - DISTRIBUTION FUNCTION FOR FREE PARTICLES 



u 




4m 2 / AS 5(p 2 - m 2 )e(p°) 8\k - p)f cd {x,p) = Am 2 8(k 2 - m 2 )e(k°)f cd (x,k) 



(65) 
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Acting similarly on (jlip we get: 

v(-k) c W-(x, k)v(-k) d = -4m 2 S(k 2 - m 2 )9(-k a )f cd (x, -k) (66) 
We now need the usual expansion of the free Dirac field (with the adopted normalizations, see sect. Hi]) : 

*(*)=E (2^37? / 5 (u r (p)e-^a r ( P )+v r (p)e^bi( P )) 
from which we can calculate the covariant Wigner function according to formula (|9|): 
W(x,k) AB = J d 4 y e- ik -y(;* B (x + y/2)* A (x-y/2) :) 

= J2WI E 7^/^< : (^b')e lp '- ( ^ /2) at( P ') + ^(p')e- ip '-^^ /2) ^') 

r,s— 1 

x f ^ (u rA (p)e- ip < x - y Ma r (p) +v rA (p)e ip <* 



2e 

Integrating the ?/ variable, the last expression becomes 



2 



(u r (p)u s (p% at(p')a r (p) :)e^'-^5 4 (k-p'/2-p/2) + v r (p)u s (p')(: a\(p')bt(p) :)c^'+^^ 4 (fc - p'/2 + p/2) 
+u r (p)i5 a (p')(: 6.0/)a r (p) :)e-^' +p )-^ 4 (fc+p72-p/2) + i; r (p)u fl (p')(: &.(p>J(p) Oe^'"^^^ + p'/2 + p/2) 

If we now integrate in d 3 x, we get: 

J d 3 *W(x,k)= j^fj ^ (u r (p)u s (p)(: ai(p)a r (p) :)<5 3 (p - p')<5 4 (fc - p) 

+ Ur (p)S s (p')(: at(p')6t(p) :)<5 3 ( P + P ')<5 3 (k + P )5(fc )e 2 ^ 0t + Ur (p)« s (p')(: & s (p'K(p) :)£ 3 (p + p')<5 3 (k - p)£(fc )e- 2 ^ 
+ Vr (p)5 s (p)(: 6 s (p)&t(p) :),5 3 (p-p')^(fc+p)) 

The two terms with (5(fc°) are contact terms which can be neglected. Integrating again in p': 

f d 3 xW(x,k)= J2 / 1 KW« S (P)(: 4(PKW :)« 4 (^P) + MpI^W ^ 4 (H)))) 

r,s — 1 

2 

= E / d 4 P^% 2 -™W)Mp)s s (p)^^ 

r.s— 1 " 

i 2 

= <5(fc 2 - m 2 )— £ [0(fc°)«r(*)u»(*O<: aJ(fe)a r (A;) :) - 9(-k°)v r (-k)v s (-k){: b\.{-k)b s {-k) :)] 

r.s—1 

Finally, multiplying by 0(±fc°) and the spinors u, v, we get: 

d 3 x u(k) c W+(x, k)u{k) d = 4m 2 9(k°)S(k 2 - m 2 )^—(: al(k)a c (k) :) 

2ek 

d 3 x v(-k) c W~(x, k)v(-k) d = -Am 2 9{-k°)8{k 2 - m 2 )^—{: b\{-k)b c {-k) :) 

2e k 

Comparing (|oT|) with (p5|) and we obtain the equation (fT2"|) . 
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APPENDIX B - CANONICAL STRESS-ENERGY TENSOR OF THE FREE DIRAC FIELD 

We want to prove eq. (fT5|) . Let us start from: 

^Vp* (tr 2 /(x, p) + tr 2 f(x, p))= [ d 4 k tr{j^W(x, k))k v (67) 



which follows from (|L3|) and (|T4|) with A = 7^. The right hand side of the above equation can be also written, by 
using the covariant Wigner function definition ([9]) as: 



J d 4 k tr(7* W(i,fc))fc" = (2^j4 E / d ^ ^T u e ~ ik ' V < : * B 



(a: + y/2)tf A (a; - y/2) :) 7lM (68) 



The integral can be worked out by parts: 



d 4 k I dZn"e- tk - v y B {x + y/2)V A (x-y/2)- I d 4 k I d 4 y e~ iky ^- (* B (x + y/2)^ A (x - y/2)) 



J d 4 k J d 4 y ^-e- ik «V B (x + y/2)* A (x - y/2) 

- f d 4 k f H 4 „p- ifc i 

dy» 

= (2tt) 4 J dS n v 8\y)^ B {x + y/2)* A {x - y/2) - (2tt) 4 J d 4 y S\y)-^- {*b(x + y/2)V A (x - y/2)) 

= - {d u ~$ B (x)* A {x) - * B (x)d v * A (x)) 

Therefore, by plugging this solution in (I68[) : 



1 , t-/ 



d 4 /c tr( 7 "W(a;, fc))fc" = ^(: *(a;) 7 ^a iy *(x) :) 



which, together with (|67p . proves eq. (TT51) . 



APPENDIX C - DISTRIBUTION FUNCTION FOR ANTIPARTICLES OF SPIN 1/2 



The goal of this Section is to prove that eq. (|22[) is in fact the Boltzmann limit of single-particle distribution function 
of antifermions with spin 1/2. The task to be accomplished is to work out the expression: 

- [V( P )eM-u/T^ z ]V{p)} T 

Henceforth, we can omit the D 1 / 2 symbol, as we work in the fundamental representation of SL(2,C). According to 
V(p) definition in eq. (fl~9|) . and using the definition of S 2 (see eq. (f2~Tj) and the equality R&(ioj/T) = exp(cj /T03/2): 

- [V(p)expl-u/TZ z ]V(p)] T = -V(pfeM-u/TZ z ]V(p) T 

- m(C \p\ ,C [p\ )^ Q R^oj/T)- 1 ) \[p]^C^ T )) 

= -m (C- 1T [p] T R (i (io;/T)- 1 [p]- ir Ct r + C T [p]^ 1T R^ /T)- 1 ^ C^ 1T ) (69) 
We can now take advantage of the properties of C = ia%, namely: 

C* T = C** = C- 1 = -C CAC- 1 = A- 1T MA e SL(2, C) 

so that the eq. becomes: 

-m (C- 1T \p} T Ra ) (iuj/T)- 1 \ P }- 1T C^ + ^[p^R^/Ty^C^) 
= m (CipfR^uj/Tr^-^C- 1 + Cb]i-^R&(iw/T)- 1 ]p]* r C- 1 ) 

= m ({pr'R^iu/Tflp} + b]tR^( lw /T) T [p]t- 1 ) = m (Ip^Rciiw/TM + [^R^/T^t" 1 ) (70) 
where, in the last equality, we have used the symmetry of the matrix R&(iu}/T). Therefore, we have shown that: 

- [V(p)expl~oj/TZ z }V( P )] T = m ([p^R^/TM + [p^R^u/T)^- 1 ) 
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whence the eq. (|22|) for antiparticles of spin 1/2 can be also written as: 

f(x,p) rs =e-te-^ i([p]" 1 R 1 ,(ia;/T)[p]) + [p]+R (a (ia;/T)[p]t- 1 )) rs (71) 

This equation is a special case (for S = 1/2) of the general expression (|18p which was obtained in ref. 14] by means 
of statistical mechanics arguments. 



